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ABSTRACT 

As planetary embryos grow, gravitational stirring of planetesimals by em- 
bryos strongly enhances random velocities of planetesimals and makes collisions 
between planetesimals destructive. The resulting fragments are ground down 
by successive collisions. Eventually the smallest fragments are removed by the 
inward drift due to gas drag. Therefore, the collisional disruption depletes the 
planetesimal disk and inhibits embryo growth. We provide analytical formu- 
lae for the final masses of planetary embryos, taking into account planetesimal 
depletion due to collisional disruption. Furthermore, we perform the statisti- 
cal simulations for embryo growth (which excellently reproduce results of direct 
iV-body simulations if disruption is neglected). These analytical formulae are 
consistent with the outcome of our statistical simulations. Our results indicate 
that the final embryo mass at several AU in the minimum-mass solar nebula can 
reach about ~ 0.1 Earth mass within 10^ years. This brings another difficulty in 
formation of gas giant planets, which requires cores with ~ 10 Earth masses for 
gas accretion. However, if the nebular disk is 10 times more massive than the 
minimum- mass solar nebula and the initial planetesimal size is larger than 100 
km, as suggested by some models of planetesimal formation, the final embryo 
mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos' 
collisional cross sections by their atmosphere could further increase their final 
mass to form gas giant planets at 5-10 AU in the solar system. 

Key Words: Planetary Formation; Planetesimals; Collisional Physics; Origin Solar 
System; Jovian planets 
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INTRODUCTION 



In the standard scenario of planetary formation, terrestrial planets and cores of Jovian 
planets are formed through the accretion of planetesimals with initial size of 10-100 km 
(e.g., Hayashi et al. 1985). This process called planetary accretion has been investigated 
by statistical simulations (e.g., Wetherill and Stewart 1993; Inaba et al. 2003; Kenyon 
and Bromley 2004), by A^-body simulations (e.g., Kokubo and Ida 1996; 2002), and by 
the hybrid method that combines an A^-body simulation for large bodies called planetary 
embryos with a statistical simulation for small bodies (e.g., Bromley and Kenyon 2006; 
Chambers 2008). Assuming that each collision of bodies leads to a perfect agglomeration, 
the growth of planetesimals can be accurately computed with A^-body codes. Inaba et al. 
(2001) showed that the evolution of mass distribution and velocity dispersion of bodies 
calculated by A^-body simulation were reproduced by the statistical simulation applying 
the collision rate and the velocity-dispersion-evolution rate based on the results of orbital 
integrations. Therefore, statistical simulations are reliable unless planetary embryos collide 
with each other after developing a long-term orbital instability. 

Planetary embryos initially form through the runaway growth (e.g., Wetherill and 
Stewart 1989; Kokubo and Ida 1996). The embryos keep their orbital separations and 
grow through collisions with surrounding planetesimals (Kokubo and Ida 1998). At this 
stage, planetesimals of almost initial size dominate the total mass of bodies (surface 
density). When embryos reach about the mass of Mars, the velocity dispersion of 
planetesimals is increased by the gravitational scattering at the embryos. Accordingly, 
another growth regime sets in, referred to as oligarchic growth. At this stage, higher 
velocities of planetesimals cause their collisional fragmentation. After a chain of successive 
destructive collisions, often called "collision cascade", bodies get smaller and smaller until 
they are removed by gas drag in protoplanetary disks or by radiation pressure and/or 
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Poynting-Robertson drag in debris disks. As a result, collision cascade decreases the surface 
density, which slows down the growth of planetary embryos (Inaba et al. 2003; Kenyon and 
Bromley 2008). 

The planetary core (embryo) exceeding the critical core mass as large as 10 Earth 
masses can no longer retain a hydrostatic envelope, resulting in the gas accretion and 
formation of the gas giant planets (e.g., Mizuno et al. 1980; Bodenheimer & Pollack 1986; 
Pollack et al. 1996; Ikoma et al. 2001). However, since the embryo growth is hindered by 
the planetesimal depletion in collision cascade, the final embryo can hardly grow beyond 
the Mars mass and thus cannot form a gas giant (Kobayashi and Tanaka 2010). 

However, it is possible that fragmentation inhibits the embryo growth to a lesser extent 
than ascertained before. Collisional fragments orbiting a central star drift inward by gas 
drag. The drift time shortens as bodies become smaller by collision cascade, until at a 
certain size their motion gets coupled with gas motion. Such coupled bodies have lower 
drift velocity and thus can survive around planetary embryos for longer time. Kenyon and 
Bromley (2009) suggested that accretion of those coupled bodies by planetary embryos may 
promote further growth of the latter. In their simulations, they assumed that the coupled 
bodies no longer experience an inward drift. However, the drift, although at a reduced rate, 
is important to determine the embryo mass gain due to the accretion of coupled bodies. 

This paper investigates the embryo growth taking into account the accretion of 
fragments resulting from the collision cascade before their removal by gas drag. To this 
end, we perform an analytic study that extends the model of Kobayashi and Tanaka (2010) 
by including the removal by gas drag, as well as statistical simulations. In the numerical 
treatment, we do not neglect the drift of coupled bodies by gas drag. Instead, we take into 
account that the gas drag law changes for small bodies in the both analytical and numerical 
procedure (e.g., Adachi et al. 1976). The goal is to find out to what extent the collisional 
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fragmentation combined with gas drag would affect tlie embryo growtli and wlietlier an 
embryo can reacli tlie critical core mass. 

We develop the analytic theory and derive the final embryo mass with fragmentation 
and gas drag in Section 2. In Section 3, we check the formulae for the final mass against 
the statistical numerical simulations. Section 4 contains a summary and a discussion of our 
results. 

2. THEORETICAL MODEL 

In this section, we summarize the oligarchic growth of planetary embryos and the 
surface density decline resulting from fragmentation. Then we derive the general formulae 
for the final mass through the oligarchic growth with fragmentation. 

2.1. DISK AND FRAGMENTATION MODEL 

We introduce a power-law disk model for the surface mass density of solids and gas. 
The solid surface mass density is taken to be 

Ss,o = /iccSi (^Y^j gcm-^ (1) 

where a is the distance from a central star, Si is the reference surface density at 1 AU, 
q is the power-law index of the radial distribution, and /ice is the factor that represents 
the increase of solid density by ice condensation beyond the snow line a^ce at which the 
temperature equals the ice condensation temperature ~ 170 K. We set the gas surface 
density to 

Sgas,0 = /gasSl (^^^j gCm"2, (2) 
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where /gas is the gas-dust ratio. In the nominal case, we adopt /gas = 240 (Hayashi 1981^ 
If the disk is optically thin, the gas temperature is given by 



where L* and Lq are luminosities of the central star and the sun, respectively. For L^ = Lq, 
this yields Oice = 2.7 AU. In the minimum-mass solar nebula (MMSN) model. Si = 7g/cm , 
q = 3/2 and /ice = 1 (a < Oice) and 4.2 (a > Oice). However, since a large amount of small 
dust is present even after completion of planetesimal formation, the disk is expected to be 
optically thick. This would make the temperature lower and its radial profile different (e.g., 
Kusaka et al. 1970). However, these effects would not drastically influence the embryo 
growth except for the location of the snow line. 

We take the fragmentation model described by Kobayashi and Tanaka (2010). 
Assuming that the fragmentation is scaled by the energy, the total ejecta mass nie produced 
by one collision between mi and m2 is given by 

^^ = ^. (4) 

rrii + 7712 1 + 

Here the scaled impact energy is given by mim2v'^ /2{mi + m2)'^Q^, where v is the 

coUisional velocity between mi and m,2 and Q^ is the critical specific impact energy needed 

to disrupt the colliding bodies and eject 50% of their mass (me = (mi + m2)/2). The value 

of Qd is given by the larger of the two colliders {m,i and m,2). Note that, since Kobayashi 

and Tanaka (2010) separately determine the ejecta mass from mi and m2 in their analysis, 

Eq. (jlj) is different from their definition. 



The energy threshold is given by 



/3= / . \/3g 2Gm 



Qh = Qos(- — +QogP(-. — +Cgg , (5) 

\lcm/ \lcm/ r 

where r and m, are the radius and mass of a body and p is its density. Benz and Ashpaug 

(1999) provide Qos, Ps,Qog, and /3g for r = 1-10^ cm from the hydrodynamical simulation of 



the collisional dispersion. The first term in the right-hand side of Eq. ([5]) controls Q^ for 
r < 10^-10^ cm and the second term describes Q^ of the larger bodies. For r > 10^ cm, 
Qd is purely determined by the gravitational binding energy, being independent of material 
properties. The collisional simulation for the gravitational aggregates yields Cgg ^ 10 
(Stewart and Leinhardt 2009). 

2.2. ISOLATION MASS 

Planetary embryos no longer grow after they have accreted all planetesimals within 
their feeding zones. The width of a feeding zone is equal to the orbital separation of the 
embryos, b{2M/3M^:Y'^a, where M is the embryo mass, M* is the mass of central star, 
6 ~ 10 is a factor (Kokubo and Ida 2002). Therefore, the maximum mass or "isolation 
mass" satisfy Miso = 27ra^(2Miso/3M*)^/^6Ss. It can be expressed as 

where M^ is the Earth mass and Mq is the solar mass. The planetary embryo mass 
approaches the isolation mass if fragmentation is ignored (Kokubo and Ida 2000; 2002). 
However, if fragmentation is included, the final embryo mass is expected to be smaller. 

2.3. PLANETARY GROWTH 

2.3.1. GROWTH WITHO UT A CCRETION OF FRA GMENTS 

At the runaway stage of planetary growth, the larger planetesimals grow faster than 
smaller ones. The solid surface density at the runaway growth stage is determined by 
relatively small planetesimals. At the later, oligarchic stage, planetary embryos become 
massive, start to gravitationally stir up the planetesimals and induce their collisional 
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disruption. The resulting fragments are quickly removed by the inward drift due to gas 
drag. Thus the fragmentation would reduce the surface density and hence the final embryo 
mass. Here, taking into account the depletion of planetesimals due to the fragmentation, we 
estimate the final embryo mass. We neglect the accretion of fragments onto the embryos, 
whereas the accretion of fragments will be considered in Section 2.3.2. 

At the oligarchic stage, planetary embryos are distributed uniformly due to their 
mutual gravitational interaction. Since they cannot collide with each other due to their 
large orbital separation, they grow slowly. If a planetary embryo with mass M collides with 
planetesimals of mass m, its growth rate is given by 

dM 



dt 



/ dmmnsa^hl,M{Pco\)^K, (7) 



where Usdm is the surface number density of planetesimals with mass in the range 
of [m,m + dm], ^k is the Kepler angular velocity, hm,M = [{'m + M)/3M^,]^^^ is the 
dimensionless reduced Hill radius between m and M, and (Pcoi) is the dimensionless 
collision rate. We assume m <^ M, and Eq. ([7]) reduces to 

—TT = C'accSga h^i{Pcol)^K, (8) 

where Eg is the surface density of planetesimals, Hm = (M/3M^,)^/'^, and C^cc is the 
correction factor on the order of unity. Because the eccentricity dispersion e* and inclination 
dispersion i* of planetesimals are much larger than Hm for kilometer-sized or larger 
planetesimals, the dimensionless collision rate (Pcoi) is given by (e.g., Inaba et al. 2001) 

(^col) = ^^, (9) 

where Ccoi — 36 for e* = 2i*. Here, R = {9M^:/Anp)^^^ /a is the embryo radius scaled by 
hMd and independent of the embryo mass. In this paper, we do not take into account 
the enhancement of (Pcoi) due to atmosphere which an embryo is expected to acquire. 
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Such an enhancement may be efficient for M > O.IM® (Inaba and Ikoma 2003). Note 
that embryos grow not only through colhsions with swarm planetesimals but also with 
neighboring embryos because their separation scaled by the Hill radius decreases by their 
growth. Furthermore, taking into account mass distribution of planetesimals, the actual 
accretion rate becomes larger. In this paper, we apply Cacc — 1-5. 

Since dM/dt oc e*~^, we need to describe evolution of e*. Two mechanisms that affect 
e* are viscous stirring and damping by gas. The viscous stirring of planetary embryos raises 
the random velocities of planetesimals. The stirring rate of e*^ is given by 

'de*^\ 
— = nMa'hUPYs)^K. (10) 

at J yg 

Since embryos are uniformly distributed with a separation 'T^l^hyiah^ their surface number 
density um is given by 

nu = — ; — • (11) 

For e* = 2i* ^ Hm, the dimensionless stirring rate (-Pvs) is given by (Ohtsuki et al. 2002) 

(Pvs) = ^^ln(A^ + l), (12) 

where Cys — 40. Although ln(A^ + 1) weakly depends on e*, we assume ln(A^ + 1) ^ 3 for 
this analysis. On the other hand, e* is damped by gas drag and the damping rate is given 
by 



dt J T 

/ gas 

where Cgas — 2.1 (e.g., Inaba et al. 2001) and 

^=—IF. • (14) 

Here, pgas is the gas density in the midplane, fK is the Keplerian velocity, and Cd is the 
dimensionless gas drag coefficient. The latter is defined as a factor that appears in the 
expression for the gas drag force acting on a planetesimal with radius r: C-QTir'^PgasU^ /2 with 
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u being the relative velocity between the planetesimal and gas. For kilometer-sized or larger 
planetesimals, Cd is a constant (Adachi et al. 1976). We assume that e* is determined by 
the equilibrium between the stirring and the damping (Eqs. fITU]) and flT5]) ). Equating the 
stirring and damping rates, we obtain 






Cvs ln(A2 + l)nKr 



nl/5 



2^/Hba 



gas 



(15) 



We can estimate e* ~ (ri^K)^^^/?.^/. Since tQk ^ 1 for the planetesimals (Adachi et al. 
1976), our assumption of e* ^ Hm for the derivation is valid. Inserting Eqs. ([9]) and flTSl) to 
Eq. ([8]), the growth rate is found to be 

1 2/5 

(16) 



dM 



dt 



(-^arr^ftCL Cfol-fl'' 



2''/^7ibC, 



Cvsln(A2 + l)fiKr 



'I'M- 



The time evolution of a planetary embryo in the oligarchic growth mode, neglecting 
fragmentation, is summarized in Appendix A. 

As embryos grow, the random velocity of planetesimals increases, making collisions 
between planetesimals destructive. Fragments produced by the collisions get smaller and 
smaller by successive collisions (collision cascade) until the smallest ones are brought inward 
by gas drag and are lost to the central star. Planetesimals with mass m dominate the 
surface density Sg during the oligarchic growth. Thus, the collision cascade induced by 
embryo growth reduces the surface density. The evolution of Sg due to a collision cascade 
is given by (see Kobayashi and Tanaka 2010; Appendix B for a derivation) 



dt 
where ho = l.lp~^/^ and 



S2(2-a)2/iofiK 



m 



1/3 



Sl23iaj 



1 + 



In 



v[m) 
2Q^(m) 



+ 



a-l 



1 



Sl23(tt), 



+ ln(l + 



(17) 



(1^ 



[1 + 4))^ ' 2-6 

where b is the power-law exponent of the mass distribution of ejecta yielded by one collision 
between mi and 777-2 and e is a factor in the expression for the mass of the largest ejecta: 
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rriL = e{mi + 1712)4'/ {1 + </')^ with a use of the scaled impact energy (p. For a coUision 
cascade, the mass distribution exponent a of fragments is given by a = (11 + 3p)/{6 + 3p) 
for f (771)^/(5^(^1) oc m^'P and S123 is insensitive to constants e < 1 and b = 1-2 (Kobayashi 
and Tanaka 2010). We set b = 5/3 and e = 0.2. 

Dividing Eq. (TTB]) by Eq. (1171) and integrating, we obtain the final embryo mass, 
^ _ /2a-l^^/(--Va^accC.oli?mV3V/^^"-V4^3a-)/(--) 



'c 



x(3M, 



3 J \ (2 - a)2/7oSi23 ; \2Qh 

Ss,0 



3/(2a-l) 
-(4-2Q)/{2a-l) 



In 



6Q:/5(2a-l) 

Cvsln(A' + l)l^Kr 



(19) 



where v = e*VK- This equation is valid unless the surface density is much smaller than that 
of embryos. Eg ^ McUm- Therefore, an estimate of Sg for the final embryo mass is given by 

cAt7^\ . (20) 



Ss,0 24/37r6a2/7M,Ss,0 AMso/ 

where Cs^ <^ 1 is a constant. Generally, Mc should be smaller than Mjso. 

We assume Qd = Qogpr'^^ for the gravity regime and thus a = [11 + 3(/3g — 2/15)]/[6 + 
3(/3 — 2/15)]. For the minimum-mass solar nebula (Eq = 7g/cm and q = 3/2), Eq. ( TT9|) 
can be re-written as 



Y ln(Ss,o/] 



/^ \ \ 1.21 / \ 0.63 / \ 0.48 

V 4.5 J \5AUJ V4.2xl020g^ 



,^ X -0.28 / ^ \ 0.89 / \ 1.85 



H^J I2.1ergcm3/g2j [Y^J^ ) '^^' ^''^ 

where /3g = 1.19 for ice. Here, we estimate ln(Ss/Ss,o) — 4.5 from Eq. ( I20l) with the use 
of Cs, = 0.1 for Mc = O.IOM® and Mjso = 2.8Me. With fragmentation, the final embryo 
mass given by Eq. (I2T1) becomes much smaller than the isolation mass, Eq. ([6]). As we 
will show in Section [3] with the aid of statistical simulations, the planetesimal mass m is 
about 100 rriQ, where mo is the initial planetesimal mass. We assume r?i = 100 rriQ and will 
compare Mc with the embryo mass obtained through the statistical simulation. 
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2. 3. 2. GROWTH WITH A CCRETION OF FRA GMENTS 

For the steady-state mass distribution of collision cascade (see Appendix A), the surface 
density of small fragments is much lower than that of planetesimals and the accretion of 
fragments is insignificant for embryo growth. The collision cascade may terminate at a 
certain mass where destructive collisions no longer occur due to low coUisional velocities 
damped by a strong gas drag. The steady-state mass distribution is achieved if fragments 
at the low-mass end of the cascade are quickly removed by the gas drag. However, when 
the Sg-decay time resulting from the collision cascade is shorter than the removal time of 
the fragments by gas drag, the termination of the collision cascade yields a large amount of 
fragments at the low-mass end of the cascade and would dominate the solid surface density. 
Here, we examine the termination of the collision cascade and estimate the embryo growth 
through accretion of such fragments. 

For small fragments, the e*-damping rate resulting from the gas drag is given by 
(Adachi et al. 1976) 

^^*^ ^*^ (22) 



dt J r T 

/ gas,f 

where rj = (vk — Vgi^)/vK is the deviation of the gas rotation velocity fgas from the Keplerian 
velocity. For q = 3/2, i] = 1.8 x 10~^(a/l AU)^'^. For small fragments, Co depends on 
u and e* becomes much smaller than rj by the gas drag. Hence, the damping rate given 
by Eq. fl22]) is different from Eq. flT^ . Since the gas drag substantially damps e* of small 
fragments, (-Pys) = 73 is independent of e* for e* <^ Hm- Thus, the equilibrium condition 
between the stirring (Eqs. (fTOl) - (fTT!) ) and the damping (Eq. (l22l) ) gives 



*2 hlj{Pys)rnK 



Smaller fragments have low e* because of low r. 

The collision cascade will no longer operate for such low e*. The collisional energy 



(23) 
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between bodies of mass m is estimated to be me*^f^/4 and is much smaller than energy 
threshold 2mQ^ for mg = m at the low-mass end of the cascade. Therefore, at the low-mass 
end 

e*'vl = C^Q*^, (24) 

where Cl is a constant. As we will show, the surface density of fragments becomes higher 
at the low-mass end of the collision cascade. Our simulation yields Cl — 1 at the low-mass 
end. 

The surface density of the fragments decreases as a result of the radial drift by gas 
drag. From Eqs. (12^ and (^^, we estimate the scaled stopping time fgtop of the fragments 
at the low-mass end of the cascade to be 



'stop 



V hlj{Pys) 4 



(25) 



^ ^^^\5A\j)[o.ImJ V3.1xl0°erg/gj' ^^^^ 

where Q^ ~ 3.1 x lO^erg/g for ice bodies with r = 10 m and we will derive M ~ O.IM^ 
in the following analysis. Since fstop ^ 1, the coupling between the fragments and gas 
is weak. Therefore, the drift velocity of fragments is given by 2rfa/T (Adachi et al. 
1976). Since r oc a?/'^ for the fragments and Eg oc a"', the Sg-depletion rate is given by 
dT,s/dt = —2(9/4 — g)r7^Ss/r. Eliminating r of fragments at the low-mass end of the cascade 
by Eq. (!25l) . we get 



dt V4 7 2V3CLvr6Qj) 

Since (Pcoi) — 11.3V -R is independent of M for e* ^ Hm, integration of Eq. (^ divided by 
Eq. fl27|) gives another formula for the final embryo mass, 

-1 3/4 



M{ 



3(9/4 - g)r/Pvs 4 ''° 



(28) 
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where we set M = at Ss = Ss,o- We recall that (-Pvs) = 73 and (Pcoi) = 11.3v i? in 
Eq. ( l28l) are independent of Mf for e* <^ Km- In addition, since QJ^ for the strength regime 
is almost independent of the fragment mass m, we neglected the mass dependence of Qd for 
the derivation of Eq. fl28l) . Assuming Ss,o ^ ^s and q = 3/2, we obtain 



3/4 / ^ ^ 3/2 / ;,^ X 3/8 

Mf = 0.14 f 



by f a Y'fM, 

loy V5Auy \Mq 
( „ y/y s,„ y/y Q° W ram 

Vlg/cmV V2.7g/cmV V3.1 x 10«erg/g/ ®' * ' 

The final mass given by Eq. (129|) is also much smaller than the isolation mass, Eq. ([6]). 
Furthermore, Mf is comparable to M^ given by Eq. fl2Tl) . When planetesimals have a large 
Qd because of a rigid material (high Qog) and/or a large size, M^ is larger than Mf. In the 
opposite case, Mf exceeds M^. The final mass is supposed to be the larger of the two, M^ 
and Mf. We apply QJJ, for r = 10 m and will compare Mf with the embryo mass resulting 
from the statistical simulation. 



3. NUMERICAL SIMULATIONS 

3.1. METHOD 

Many authors attacked the problem of the planetary growth from planetesimals to 
embryos by A^-body simulations (e.g., Kokubo and Ida 1996; 2002), the statistical method 
(Wetherill and Stewart 1993; Kenyon and Bromley 2004), and the hybrid method (Kenyon 
and Bromley 2008, Chambers 2008). We apply the statistical method developed by Inaba 
et al. (1999, 2001), which we briefly explain now. The mass distribution of bodies in orbit 
around a central star evolves through mutual collisions, taking into account gravitational 
focusing of colliding bodies (see Inaba et al. 2001 for detailed expressions). Along with 
the mass evolution, velocity dispersion changes by gravitational perturbations, coUisional 
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damping, and gas drag. Equations for the mass distribution of bodies are integrated 
simultaneously with equations for their velocity distribution. 

The time evolution of the differential surface number density ns{m) at a distance a is 
given by 

dmns{m,a) m. 



oo 



= — JIk / dmi / dm2 

^ ^ Jo Jm—mi-rnc 

■x{hra,,m2a)^'n's{mi,a)n^{m2,a){Pcoi) 

x6{m — nil — '"^2 + '^e) 

-nKmrisim) / dm2{hra,m^afns{m2,a){Pco\) 
Jo 

Q POD pmi 

+ -— r^K / dmi / dm2{mi + 1712) f{m, mi, 1712) 

xrisinii, a)ns(m2, a)(/imi,m2a)^(Pcoi) 

1 d 

— ■—[amns{m,a)vdriit{m,a)], (30) 

a oa 

where S{x) is the delta function, {mi + m2)f{m,mi,m2) is the mass of fragments less than 
m produced by a collision between irti and m2, and irte is the total mass of the fragments 
(given by Eq. (jl]) in our fragmentation model). We apply the dimensionless coUisional rate 
(Pcoi) which Inaba et al. (2001) provide as a function of mi, 777-2, and e* and i* by compiling 
previous studies. Equation (!30|) describes the mass transport in the two-dimensional 
space composed of mass and distance (radial direction). The first and second terms in 
the right-hand side of Eq. (130 p represent the mass transport along the mass coordinate 
caused by coagulation and the third term does that due to fragmentation. The fourth term 
describes the mass transport along the radial coordinate due to the drift of bodies. We 
calculate the transport on a grid of mass and radial bins. 

Assuming a power-law mass distribution of fragments, (7771 + 7772)7(777,7771,7^72) is given 



by 
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■rric I 1 for m < m-L, 

{mi+m2)f{m,mi,m2) = { \^lJ (31) 



j5 

\2 



rUe for m > rrii^, 

where we recall the mass of the largest ejecta rriL = e(mi + m2)0/(l + 0)^ and the total 
ejecta mass rrie given by Eq. @. 

The drift velocity is characterized by the dimensionless parameter fstop, where the 
scaled stopping time fstop is given by aQ'^r/u. The relative velocity u between a body and 
gas is equal to (e* + i* + ri)aQK- For fgtop ^ 1 and fgtop ^ 1, the drift velocity can be 
written as (Adachi et al. 1976, Inaba et al. 2001) 

'"""■^ '-°- for f..„„<l, 



~2 '■^^ 'stop r^ ^1 

Vdnlt[m,a)={ ''*°''. r„^j'.„ ,.. _.„,.xnO , .1/2 



av r mm+Kimi ^ 1.2 ^ 2 ^ '^" ^^^^ 



''"stop -L + ''"stop 

where E and -ft' are the elliptic integrals. Both regimes in Eq. ( I32l) can be combined into 

''"-'- + ''"stop L Jvr vr J 

Indeed, if fstop ^ 1, then e* and i* are damped by the strong gas drag to e* ^ r^ and i* ^ r], 
and Eq. fl55]) reduces to the first of Eqs. fl52]) . If fstop ^ 1, Eq. fl5^ simply coincides with 
the second of Eqs. ( 132|1 . 

Since we treat small fragments, we consider three gas-drag regimes: quadratic, Stokes, 
and Epstein ones. The gas drag coefficients Cd for these regimes are given by (Adachi et 
al. 1976) 

Cn.quad = 0.5, Cd, Stokes = 24:/Re, Cd, Epstein = l6/KRe, (34) 

where K = 1.7 (10~^gcm~^/pgasr) is the Knudsen number and Re = AA{u/c)K^^ is the 
Reynolds number with c being the sound velocity of gas (Adachi et al. 1976). All three 
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regimes can be described together by adopting 

Cd = C'o.quad + I -^ h 7^ 1 • (35) 

\ L^D, Stokes <-^D,Epstein / 

Since the mean colhsion rate is a function of e* and i*, the evolution of e* and i* 
should be calculated precisely. We consider gravitational perturbations from other bodies, 
collisional damping, and gas drag. Their net effect can be calculated as a square of 
dispersions: 

+ -IT ' (36) 



^^ V c?t /g,,, V c^t ) ^^ \ dt /^„u 



where terms with subscripts "grav", "gas", and "coll" indicate the time variation due to 
gravitational perturbations, gas drag, and collisions, respectively. These terms are provided 
by Inaba et al. (2001). Note that for the collisional damping, both fragments and mergers 
resulting from a collision have the velocity dispersion at the gravity center of colliding 
bodies, as described in their Eq. (30). We follow the planetary growth by simultaneously 
integrating Eqs. (|30|), (|36l), and (|371). 

Kokubo and Ida (1998) and Weidenschilling et al. (1997) showed that the radial 
separation of orbits of runaway bodies formed in a swarm of planetesimals is about 10 
times their mutual Hill radius. Since dynamical friction from the field planetesimals makes 
orbits of the runaway bodies nearly circular and coplanar, the orbital crossing of these 
runaway bodies never occurs before the gas is dispersed (Iwasaki et al. 2001). Therefore, 
they cannot collide with each other. In our simulation, when the bodies reach a certain 
mass mrun such that the sum of their mutual Hill radii equals the radial-bin width divided 
by 6, the bodies are regarded as "runaway bodies" which do not have any collisions and 
dynamical interaction due to close encounters among them, following Inaba et al. (2001). 
We set h = 10. As the bodies grow, their separation measured in the Hill radii decreases 
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and these bodies can no longer be treated as "runaway bodies" . Therefore, 171^-^^ increases 
during the simulation and then the number of the "runaway bodies" decreases. 

Moreover, Inaba et al. (2001) set the number of collisions during a numerical time 
interval to be an integer using a random number generator and hence keep the number of 
bodies an integer. Our procedure is different. Instead of dealing with the number of bodies, 
we treat the mathematical expectation N of the number of bodies with mass larger than m. 
When N is much lower than unity, the bodies may not yet exist. We introduce a certain 
critical number A^c (a "threshold") to get rid of such bodies. If N{mc) = Nc, we neglect 
collisions with bodies larger than rric and dynamical interactions with them. The value of 
Nc should be of the order of unity, because we can say that bodies are not yet born for 
A^ <^ 1, whereas choosing a high N^ would delay the embryo growth. Bodies with masses 
ranging from mj-^n to rric are treated as "runaway bodies" . The mass range of the "runaway 
bodies" tends to extend during the simulation for low N^, although being kept small in 
A^-body simulation (Kokubo and Ida 2002). We always start the simulation with N^ = 0.1 
and A'c may change from 0.1 to 10 during the simulation to keep the small mass range of 
the "runaway bodies" . As shown in Figs. [1] and [21 our simulation with such a choice of 
A"c reproduces the mass and velocity-dispersion distributions obtained from the A^-body 
simulation, presented in Figs. 4 and 5 of Inaba et al. (2001). Where about 100 runs of the 
statistical code by Inaba et al. (2001) are required to reproduce the A^-body simulation, our 
simulation does the same with only one run. 
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3.2. EMBRYO GROWTH 

3.2.1. WITHO UT FRA OMENTA TION 

We calculate the evolution of the number and velocity dispersion of bodies by summing 
up the time variations coming from all mass and radial bins and by integrating them over 
time. We integrate Eq. f l5U]) through the fourth order of Runge-Kutta method for the mass 
evolution and Eqs. ( 136|) and (!37|) through the linear method for e* and i* evolution. We set 
the mass radio between adjacent mass bins to 1.2, which we found sufficient to reproduce 
A^-body simulations (see Figs. [Hand [2]). Fig. 12] shows the results ignoring fragmentation 
{Qd = ^^) with a set of six concentric annuli at 3.2, 4.5, 6.4, 9.0, 13, and 18 AU for 
M^, = Mq, Si = 7g/cm , and q = 3/2. We consider planetesimals with initial mass of 
mo = 4.2 X 10^® g (radius of Tq = 10 km for p = 1 g/cm ) and velocity dispersion given 
by e* = 2i* = (2mo/M*)^/^. The planetesimals are assumed to be composed of ice whose 
physical parameters are listed in Table 1. We artificially apply the gas surface density 
evolution in the form Sgas = Sgas,oexp(— t/Tgas,dep), where Tgas,dep is the gas depletion 
timescalqj. Here, we set Tgas,dcp = 10^ years. 

Figure [3^ shows the mass distribution of bodies at 3.2 AU from 10^ years to 4 x 10^ 
years. After 10^ years, the mass distribution of small bodies {m < 10^^ g) is a single power 
law, which is consistent with that of the runaway growth, n^ oc m^^^^ (Makino et al. 1998). 
At the same time, the oligarchic growth results in the fiatter distribution for large bodies 
{m > 10^^ g). Such large bodies have almost a single mass (~ 3 x 10^^ g at 10^ years and 
~ 2 X 10^^ g at 4 X 10^ years). As the growth of the large bodies proceeds, the random 
velocities of small bodies exceed the escape velocity because of the viscous stirring by large 



^Assuming a constant Sgas gives almost the same results for the final embryo mass, 
because we consider time spans t < Tgas,dep- 
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bodies in the entire mass range of small bodies [m < 10^^ g). Then, the mass distribution 
starts to depart from this power law. Collisions between small bodies at such velocities 
lead to a fiat mass distribution for small bodies (m < 10^^ g at 10^ years and m < 10^^ g at 
4 X 10^ years) and a steep one for large bodies. This means that planetesimals dominating 
the surface density become larger after 10^ years. 

The velocity dispersion evolution (e* and i*) at 3.2 AU is shown in Fig. HI At 10^ 
years, e* and i* are proportional to m~^/^, being determined by the dynamical friction. At 
10^ years, e* and i* for m < 10^^ g are independent of mass, whereas e* and i* of larger 
bodies still preserve an m"^^'^ dependence. The viscous stirring by large bodies dominates 
e* and i* of small bodies. In this case, e* is given by Eq. ( 1A1|) . After 10^ years, the velocity 
dispersion of small bodies {m < 10^^ g at 10^ years and m < 10^^ g at 4 x 10^ years) reaches 
the equilibrium between the stirring and the gas drag. In this case, e* is given by Eq. (ITSl) . 
However, bodies with intermediate mass do not yet have the velocity dispersion in the 
equilibrium. The velocity dispersion of bodies dominating the surface density is still given 
by Eq. ( JAll) . The growth of large bodies are described by M^ given by Eq. ( JASP in 10^-10^ 
years (see Appendix A). 

In the outer region, the mass distribution and velocity dispersion evolve in a similar 
way but in a longer timescale. At 4 x 10^ years, the bodies with m > 10^^ g grow through 
the oligarchic mode at a < 9 AU. Beyond 9 AU, the largest bodies have not yet reached 
10^^ g by that time. 

3. 2. 2. WITH FRA OMENTA TION 

We now take into account fragmentation, using Q^ in Eq. ([5]). Kenyon and Bromley 
(2009) suggested that the small bodies coupled with gas (< 1 m) affect the embryo growth. 
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To treat the coupled bodies, we set the smallest mass to be 4.2 g (radius of 1cm for 
p = Ig/cm ). 

The mass distribution and velocity evolution that we calculated with fragmentation 
are shown in Figs. [5] and [61 During the first 10^ years, bodies grow through mutual 
collisions. Fragments {m < mo) are not yet numerous because of the low velocity dispersion 
(Fig. [6]). At t = 10^ years, the runaway growth occurs inside 6.4 AU. Similar to the case 
without fragmentation, the cumulative number Uc is almost proportional to m~^^^ for 
m > 10^° ~ IOOttiq and the distribution of large bodies with m > 10^'^ g is flat. Since the 
velocity dispersion of planetesimals ~ e*f k exceeds their escape velocity for m < 10^"^ g, 
planetesimals dominating the surface density become larger after 10^ years. In contrast to 
the case without fragmentation, the mass of bodies dominating the surface density ceases 
at about lOOrrio, owing to a high collisional velocity e*VK ^ \/Qb- The slope of ric for 
m ~ 10^^-10^'^ g ('^ = 0.1-10 km) is nearly —(5 + 3p)/(6 + 3p), which is typical of a collision 
cascade for v^/Q^ oc m~P (Kobayashi and Tanaka 2010). The downward mass flux along 
the mass coordinate is constant with mass in the collision cascade. For m < 10^^ g, the 
velocity dispersions (e* and i*) are effectively damped by the gas drag in the Stokes regime 
(See Fig. [6]). Because the collisional energy (~ e*^f^) is lower than the energy threshold 
(Qd) for a destructive collision, the downward mass flux becomes insignificant compared 
to that in the collision cascade. Therefore, the number of the bodies increases around 
10^^ g where e*'^v^ ~ Q^. Bodies with m ~ 10^ g (r ~ Im) are most efficiently removed 
by the radial drift and smaller bodies {m ^ 10^ g or r ^ 1 m) are coupled with gas. As a 
result, the mass distribution becomes wavy (see Fig. |5]). The surface density of bodies with 
e*'^v^ ~ Qd is much higher than that of the coupled bodies. However, planetesimals with 
~ lOOrrio dominate the surface density. At 4 x 10^ years, the oligarchic growth is inhibited 
by fragmentation inside 5 AU, where the number of fragments decreases compared to that 
at 10^ years. 
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Figure [7] shows the embryo-mascl evolution at 3.2 AU. At the runaway growth phase 
(the first < 10^ years), the embryo mass grows as an exponential function of time. The 
growth rate becomes much lower in the oligarchic mode at 10^-10^ years. The embryos 
exceed that for the case without fragmentation. They are substantially fed by small 
fragments whose velocities are damped by gas drag and are low (see Figs. [6]), resulting in 
their rapid growth. However, they grow at a sluggish pace at about 4 x 10^ years because 
the surface density of planetesimals is decreased by collision cascade followed by removal of 
small bodies by gas drag. As a result, the final embryo mass is reduced by fragmentation. 
In this case, since M^ and Mf are comparable, the embryo mass is determined by both 
effects. Even so, the final embryo mass is consistent with M^ and Mf. 

Fig. [8] presents the embryo mass at 10^ years as a function of distance, including the 
results for three additional annuli at 1, 1.4, and 2.0 AUo In simulations at a < 2.7 AU, 
we use silicate bodies (Table 1) instead of icy bodies. Decline of the planetesimal surface 
density caused by fragmentation results in embryos smaller than those for the case without 
fragmentation inside 6AU. The masses reach only 0.01-0. IM^ in 10^ years. The embryo 
masses at a < 1.5 AU and a = 2.7-5 AU are consistent with Mc or Mr. Inside a < 1.5 AU, 
the final mass is determined by the largest of the two masses, Mc and Mf (see Fig. [9]). 
For rriQ = 4.2 x 10^^ g (rg = 1km), the embryo mass is determined by Mf because of 
Mf > Mc (see Fig. [10]). In addition, if both M^ and Mf are smaller than M^ in Eq. flA3l) . 
the fragmentation is nearly negligible and the embryo mass is given by Mn (see Appendix 
A for an analysis of embryo growth without fragmentation). 

Since the initial mass of planetesimals depends on their formation process, which is not 
well understood, it is worthwhile investigating the embryo mass dependence on the initial 



^Here, embryo mass represents the average mass of the "runaway bodies" . 

■^The results come from two simulations: one in a < 2.7 AU and the other in a > 2.7AU. 
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planetesimal mass. Fig. [TT] shows the dependence of the embryo masses on rriQ at 3.2 AU 
for the minimum mass solar nebula after 10 million years. If the initial mass is smaller 
than ~ 10^^ g, the embryo mass is determined by accretion of small fragments. For larger 
mo, growth is halted by the collision cascade. The embryo mass increases with mo because 
the timescale of the surface density decline due to collision cascade is longer for larger 
planetesimals, strengthened by self-gravity (large Q^). If m,o > 10^^ g, the embryo grows 
without substantial fragmentation. For a massive disk (Si = 70g/cm ), the embryo mass 
can reach the Earth mass. In this case, Mc and Mf determine the embryo mass because of 
the short growth time for Tq = 1-100 km. 

4. SUMMARY AND DISCUSSION 

In this paper, we studied analytically and numerically the growth of planetary embryos 
in the framework of the standard scenario. We took into account that embryos growing 
in the oligarchic mode pump up relative velocities of planetesimals, which causes their 
collisional fragmentation. We also considered the fate of smallest fragments of the resulting 
collision cascade, namely their inward drag in the ambient gas and possible accretion by 
nascent embryos. 

Taking into account fragmentation, we have analytically derived the final mass of a 
planetary embryo in Section |5] in two cases: the mass M^ of an embryo growing through the 
accretion of planetesimals which are removed by collision cascade and the mass Mf in the 
case where an embryo grows through collisions with fragments which are removed by gas 
drag. In Sec. El we showed that the final embryo mass obtained in numerical simulations 
can be reproduced by the larger of the two analytic estimates, MAX(Mc, Mf). If the embryo 
mass Mn of Eq. flA3|) for the case without fragmentation is smaller than MAX(Mc,Mf), 
fragmentation is negligible and the final mass is given by M^. Altogether, our analytical 
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formulae provide a good estimate of the final mass of planetary embryos. 

The gas giant planet formation via core accretion requires the solid core as large as 
~ lOM^ {e.g., Ikoma et al. 2001). Kenyon and Bromley (2009) suggested that the planetary 
embryo may further grow by the accretion of bodies coupled with gas because the coupled 
bodies are no longer involved in collision cascade. Taking into account three regimes of 
gas drag (see Eq. [35]), we found that the collision cascade halts at bodies having t>^ ~ Q^, 
i.e. at sizes that are larger than the critical coupling size corresponding to fstop = 1- As a 
result, coupled bodies should be produced only in little amounts and can hardly contribute 
to planetary growth, despite their low velocity dispersion in the laminar disks considered 
here. 

In this paper, we derive the final embryo masses, assuming a uniform distribution 
of solid bodies. However, the planetary embryos can open gaps in a solid disk, which 
would affect the accretion of fragments on the embryo and their removal by the gas drag. 
Tanaka and Ida (1999) derived a condition for the gap opening around a migrating embryo. 
Replacing the embryo's migration time by the drift time T/2iif for fragments, we can obtain 
a criterion for an embryo with mass M to open a gap in the solid disk from their Eq. (3.6), 



rfix/iif 



>0.81 



1 + 0.45(rfiK/27r)2/3 + 1 



(38) 



47]-7^2 

where r is given by Eq. (TT4l) for fragments surrounding the embryos. Since Eq. (138|) requires 
tVL]^/2ti ^ 1 for M > O.IM®, the critical fragment radius fg above which Eq. (!38l) is 
satisfied is given by 

where we adopt q = 3/2 and Cb = Cc,stokes- Levison et al. (2010) performed A^-body 
simulations involving fragments and found rg ~ 30 m for M > M®, which is consistent 
with Eq. ( 139|) . On the other hand, as we have shown, the solid surface density is reduced 
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by the elimination of fragments at the low-mass end of collision cascade due to gas drag. 
In addition, the embryo grows through the accretion of the fragments for small initial 
planetesimals. The radius rg of the fragments is estimated by Eq. fl2S]) as 



re ~ 5.0 —— m. (40) 

Eqs. ( 15^ and (HO]) give r^ ^ rg. This means that the collision cascade quickly grinds 
planetesimals down to the size much less than r^. As a result, the gap formation in the 
solid disk composed of fragments with the radius r^ does not occur and the embryo growth 
would not be influenced. 

Considering the range of initial planetesimal radii tq = 1-100 km, and the range of 
disk surface densities Si ~ 7-70 g/cm at 3.2AU (1-lOxMMSN), our results suggest that 
the fragmentation averts the planetary growth at M ~ O.I-IOM0 at several AU. This 
is consistent with the results of planetary growth simulations beyond 30 AU by Kenyon 
and Bromley (2008). The cores approaching the critical mass of 10 M^ can only form at 
a = 3-4 AU, only from planetesimals that have r^ ~ 100 km initially, and only in densest 
nebula with Si ~ 70 g/cm . Thus it appears problematic to explain formation of giant 
planets at 5-10 AU, such as Jupiter and Saturn in the solar system, even for Si ~ 70 g/cm 
and Tq ~ 100 km. 

An effect that could help further is the enhancement of the embryo's coUisional cross 
section due to a tenuous atmosphere acquired by the embryo in the ambient gas. Inaba 
and Ikoma (2003) and Inaba et al. (2003) have shown that a core with mass exceeding 
about O.IMq) could grow by that effect to the critical core mass. As found here, the embryo 
masses can indeed be about O.IM^ under the minimum mass solar nebula conditions, at 
least for an initial planetesimal radius greater than 10 km. 

On any account, our results strongly suggest that increasing the disk surface density 
and/or the initial embryo size helps forming larger cores. This brings up the question 
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whether the values that we need to grow the embryo to the critical mass, lOxMMSN and 
Tq ~ 100 km — or less, if the atmospheric effect is taken into account — appear plausible. 
As far as the disk density is concerned, the answer is probably yes. lOxMMSN is close to 
the gravitational stability limit, and all densities below this limit cannot be ruled out. For 
instance, Desch (2007) and Crida (2009) point out that the MMSN profile is inconsistent 
with the "Nice model" (Gomes et al. 2004) and should be replaced with another surface 
density profile, which would possibly imply larger surface densities at several AU. 

We now consider the initial size of planetesimals. The mechanisms of planetesimal 
formation are highly debated but, despite intensive effort, remain fairly unknown. Classical 
models in which dust smoothly grows to planetesimals with ro ~ 1 km face the "meter 
barrier" problems: first, meter-sized objects should be lost to the central star as a result of 
gas drag (Weidenschilling 1977, Brauer et al. 2008), and second, further agglomeration of 
meter-sized objects upon collision is problematic (Blum and Wurm 2008). Accordingly, in 
the last years, competing scenarios were suggested that circumvent these barriers, such that 
the "primary accretion" mechanism proposed by Cuzzi et al. (2008) and "graviturbulent" 
formation triggered by transient zones of high pressure (Johansen et al. 2006) or by 
streaming instabilities (Johansen et al. 2007). These models all imply rapid formation of 
rather large planetesimals in the ro ~ 100-1000 km range. Support for these scenarios 
may come from the analysis of left-over planetesimals in the solar system. Morbidelli et al. 
(2009), for instance, suggest that the initial planetesimals should be larger than 100 km to 
reproduce the mass distribution of asteroids in the main belt. 

Our last remark is related to the Type I migration of bodies due to interaction with gas 
(e.g., Tanaka et al. 2002). The planetary embryo grows in the runaway mode followed by 
the oligarchic one. The timescale of the runaway growth is proportional to the initial size of 
bodies. In addition, the embryos rapidly grow through collisions with small fragments for 
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small Tf) even in the the ohgarchic growth. The embryo eventually grows to a larger mass if 
it is large initially, while it forms earlier if it starts with a smaller mass. Therefore, initially 
small planetesimals seem to be preferred for the core formation prior to their removal due 
to Type I migration. However, this is only valid for laminar disks. Turbulence may help 
prevent the loss of growing embryos to Type I migration (e.g., Laughlin et al. 2004). So, 
the alternative scanarios of planetesimal formation that all require turbulent disks may help 
here, too. 
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A. PLANETARY GROWTH WITHOUT FRAGMENTATION 



When the growth timescale Tgrow of embryos is much longer than the e*-damping 
timescale Tg^s by gas drag, e* settles to the equilibrium between the stirring and the gas 
drag and is given by Eq. flT^ . However, the growth timescale at the beginning of oligarchic 
growth is relatively short. In this case, e* cannot reach the equilibrium for Tgrow ^ Tgas- 
Using these two limits of e*, we will derive the embryo mass evolution as follows. 

In the non-equilibrium case (Tgrow <^ ^gas), e* increases as M grows. Then, combining 
the growth rate (Eqs. ([H])-®) with the stirring rate (Eq. (ITO|) with Eqs. (ITTl) and (IT2l)). we 
get 

(Al) 



3Cvshi(A2 + l)M4/3 



.2W3,rCacc&a2Ss/?Ccoi(3M,)i/3j 
where we assume that the initial embryo mass and eccentricity dispersion are much smaller 
than M and e*, respectively. Inserting Eq. flAip into Eq. ([8]) with Eq. (|9]) there results 



dM 



2io/37r6(CaccSsCeoia'i?)'fiK 



(A2) 



dt 9M,Cvsln(A2 + l) 

If the embryo grows with this rate for sufficiently long time, it is valid use M = at t = 
as the initial condition. In this case, the embryo mass M^ at time t is given by 



M. 



~\ 2 



2^37^6 CaccSsCcoia^i? n^t 



9M*Cvs ln(A2 + 1) 



(A3) 



which is independent of m. 



In the equilibrium case (Tgrow ^ ^gas), the growth rate of M is given by Eq. flT6|) . Also 
we consider that the embryo mass becomes much larger than the initial mass. Integrating 
Eq. flT^ over time with M = at t = 0, we have 



M. 



CaccSa Ccol-R^K 

3(3M,)2/3 



24/37r6a 
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Cvsln(A2 + l)fiKr. 
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(A4) 
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The embryo mass is proportional to m?^^ from the m dependence of r. The growth time of 
embryo mass Mg is similar to that estimated by Kokubo and Ida (2002). 

Fig. [13] shows the embryo growth at 3.2 AU. The runaway growth swithes to oligarchic 
growth after 10^ years. The embryo mass is consistent with M^ given by Eq. (1A3|) in 10^-10^ 
years. On contrary, the embryo mass does not agree with Mg given by Eq. (lA4p . That is 
because e* of the representative planetesimal, which dominates the surface density, does 
not reach the equilibrium. As seen in Fig. [3l the representative planetesimal grows because 
its random velocity is higher than its surface escape velocity. Although e* reaches the 
equilibrium for small bodies, that of the representative planetesimal cannot do due to the 
growth. Since the representative planetesimal growth is caused by the perfect sticking even 
for a high velocity, the growth came from the artificial treatment. In practice, the growth 
halts at ~ 100 mo when fragmentation is included, because the impact energy exceeds Qd- 
In this case, e* of the representative planetesimals reach the equilibrium and the embryo 
mass is determined by M^ or Mf because of the active fragmentation. 

B. MASS DEPLETION DUE TO COLLISION CASCADE 

Here, we derive the surface- density depletion rate due to collision cascade, following 
Kobayashi and Tanaka (2010). Focusing on the mass transport by fragmentation, we 
neglect terms in the right-hand side of Eq. fl5U]) except for the third one: 



dmn^ d 



poo pT^l 

VIk / dnii / dm2{Tni + 1712) /{m, mi, 1712) 

Jm Jo 

xns(mi, a)ns(m2, a)(/imi,m2«)^(-Pcoi)mi,m2 (Bl) 



dt dm ./„ ./Q 



Collisional velocities exceed the surface escape velocity of colliding bodies in collision 
cascade, resulting in (/imi,m2'j^)^(-Pcoi) = (17.3/27r)(ri + r2Y = h^m.-^ [1 + (m2/mi)^/^]^. 
Here, ri and r2 are the radii of mi and m2, respectively. Since collisions with mi ^ m2 
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mainly lead to the mass transport along the mass coordinate in collision cascade, the / 
function given by Eq. fl^ becomes 

f{m, mi, 1712 



m 
nil 



1 + 



'^ + ^f 



for ^ < jS^, 



(B2) 



1 + 



for ^ > j-^. 



Furthermore, (/^„„^,a)2(P,„,) = homf and = ivimif/2QUmi))im2/mi). Introducing 
the dimensionless variable x = m/rrii, we can re- write Eq. (jBip as 
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e4,/[l+4>? 
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1 + 



6-2 



(B3) 



where ng = Am~'^ and v'^{m)/Q'^{m) oc in~^. The upper integration limit over in Eq. (IB 3 1) 
is different from Eq. (31) of Kobayashi and Tanaka (2010), which comes from the different 
definition of /. Despite that, Eq. f lT7|l is consistent with Kobayashi and Tanaka because of 
a insignificant contribution of collisions with > v'^{mi)/2Q'^{mi) to the integral over in 
their Eq. (31). Since t>^(mi)/2Qi)("^i) ^ l, we change the upper limit to oo. We also take 
into account that v'^{mj)/Q'^{m) oc m~^ and that dm,ns/dt = in a steady-state collision 
cascade. Then, Eq. (|B3]) gives a = (11 + 3p)/{6 + 3p). Note that Eqs.(31), (34), (B2), and 
(B3) of Kobayashi and Tanaka (2010) contain mistakes: their power-law exponents of m, 
should be 11/3 - 2a, as in Eq. flB3|) . 

In the oligarchic growth, since planetary embryos sufficiently increase the collisional 
velocity of the swarm planetesimals, colhsion cascade depletes the swarm. Smaller bodies 
quickly reach the steady state of collision cascade. Therefore, integrating Eq. (1B3P over 
mass in the steady-state collision cascade (a = (11 + 3p)/(6 + Sp)), we obtain 

2tt/j-)(m,jnaxj / Jo 



dt 
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3-2 .l 

+ dxx 



-1 



(B4) 



where mmax is the mass of largest planetesimals in colhsion cascade. Integrating over x 
on the right-hand side of Eq. flB4p . and using Eg = y4m^°/(2 — a), we obtain Eq. flTTl) . 
In Eq. (TT7|) . we convert rrimax to m, according to the definition in Section [2l Moreover, 
Eq. (TT71) includes the additional terms due to the mass transport of the remnant according 
to Kobayashi and Tanaka, although the transport terms are much smaller than othersc 



''The negative sign in the right-hand side of Eq. (A5) of Kobayashi and Tanaka (2010) is 
in error, and we correct this in Eq. ([1] 
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Table 1: Material properties 
Qos Ps 

(erg/g) 



Qog 
(ergcm'^/g^) 



13, C, 



gg 



(g/cm^ 



ice 
silicate 



7.0 X 10^ -0.45 
3.5 X 10^ -0.38 



2.1 
0.3 



1.19 9 
1.36 9 



1 
3 
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Fig. [1] — Evolution of the mass distribution of bodies with mo = 2 x 10^^ g, p = 2 g/cm , 
and Cq = 2zq = 0.002 at 1 AU. The initial mass distribution is shown with the dotted curve 
in panel (a). Collisional fragmentation is neglected. Gray dashed lines represent the results 
of A^-body simulation shown in Fig. 4 of Inaba et al. (2001). 

Fig. [2] — Evolution of the velocity dispersion of bodies for the same conditions as in 
Fig. [TJ The solid and dotted lines indicate e* and i*, respectively. Circles and triangles 
represent e* and i* calculated by A^-body simulation shown in Fig. 4 of Inaba et al. (2001), 
respectively. 

Fig. [3]— The mass distribution of bodies at lO^(A), 10^ (B), 10^ (C), 4 x 10^(0) years, 
neglecting fragmentation. Different panels correspond to different radial annuli. 

Fig. H— The velocity-dispersion distribution at 3.2 AU after 10"^, 10^ 10^ 4 x 10^ 
years of evolution. Fragmentation is neglected. The solid and dotted lines indicate e* and 
i*, respectively. 

Fig. [5] — Same as Fig. [3l but with fragmentation. 

Fig. M — Same as Fig. HJ but with fragmentation. 

Fig. [3 — Evolution of the embryo mass at 3.2 AU with fragmentation (black solid line) 
and without fragmentation (gray solid line). Dashed lines indicate Mc (Eq. [191 black) and 
Mi (Eq.EH gray). 

Fig. M — The embryo mass at 10^ years for initial planetesimal mass of 4.2 x 10^^ g 
with fragmentation (filled circles) and without fragmentation (open squares). Dashed lines 
indicate M^ (Eq. [T^ black) and Mf (Eq. [2D1 grey). Gray solid line shows Miso (Eq. E]) for 
comparison. The vertical line represents the snow line a = Oice = 2.7 AU. 

Fig. [9]— Same as Fig. El but at 1 AU. 
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Fig. [To] — Same as Fig. [71 but for mo = 4.2 x 10^^ g (ro = 1km). 

Fig. [m — The embryo mass at 3.2 AU at 10^ years for the minimum-mass solar nebula 
model (El = 7g/cm and q = 3/2) as a function of the initial planetesimal mass. The 
symbols and lines are same as Fig. [81 

Fig. [121 — Same as Fig. [TTl but for Si = 70g/cm . 

Fig. [13l — Evolution of the embryo mass at 3.2 AU, neglecting fragmentation. The 
dashed lines indicate M^ (Eq. IA3t black) and M^ (Eq. IA4j gray). 
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